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Abstract 

Fully differential next-to-next-to-leading order calculations require a method to can- 
cel infrared singularities. In a previous publication, I discussed the general setup 
for the subtraction method at NNLO. In this paper I give all subtraction terms for 
electron-positron annihilation associated with one-loop amplitudes with one unre- 
solved parton. These subtraction terms are integrated within dimensional regular- 
ization over the unresolved one-particle phase space. The results can be used with 
all variants of dimensional regularization (conventional dimensional regularization, 
the 't Hooft-Veltman scheme and the four-dimensional scheme). 
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1 Introduction 



Higher order computations are a prerequisite for precise predictions in perturbative QCD. 
For many processes next-to-next-to-leading order (NNLO) calculations are required. This 
is in particular true for processes involving hadronic jets at high energies. At the next- 
to-next-to-leading order level the ingredients for the third order term in the perturbative 
expansion for quantities depending on n resolved "hard" partons are the n-parton two-loop 
amplitudes, the (n+l)-parton one-loop amplitudes and the (n+2)-parton Born amplitudes. 
Taken separately, each one of these contributions is divergent. The two-loop amplitudes 
have explicit poles up to 1/e 4 , when calculated within dimensional regularization. As usual 
the dimensional regularization parameter is e = 2 — D/2. The one-loop amplitudes have 
explicit poles up to 1/e 2 . In addition, integration over the phase space region where two 
partons are "close" to each other, brings another two powers of 1/e. Finally, the Born 
amplitudes with (n + 2) partons are finite, but integration over the phase space regions, 
where the (n + 2)-parton configuration degenerates to a (n + 1)- or n-parton configuration 
yields poles in e up to 1/e 4 . 

Due to the large variety of interesting jet observables it is desirable not to perform 
a NNLO calculation for a specific observable, but to set up a computer program, which 
yields predictions for any infra-red safe observable relevant to the process under consider- 
ation. This requires to work with fully differential cross-sections and to use Monte Carlo 
techniques for the phase space integration. 

The divergences from the separate terms discussed above cancel each other and the 
sum of the three contributions is finite. However, the cancellation occurs only after the 
integration over the unresolved phase space has been performed and prevents thus a naive 
Monte Carlo approach for a fully exclusive calculation. It is therefore necessary to cancel 
first analytically all infrared divergences and to use Monte Carlo methods only after this 
step has been performed. 

Due to the universality of the singular behaviour of QCD amplitudes in soft and 
collinear limits, the problem can be solved in a process-independent way. Possible ap- 
proaches are the extension of the phase space slicing method or the subtraction method 
to NNLO. Both methods are well understood in NLO calculations. A third and different 
approach at NNLO does not try to solve the problem in a process-independent way, but 
uses for specific observables the optical theorem to relate phase space integrals to loop 
integrals [1, 2]. The phase space slicing method has been used for the NLO calculation 
of the photon + 1 jet rate in electron-positron annihilation [3]. This calculation shares 
already some features with a full NNLO calculation. 

Here, I focus on the subtraction method. This method consists in adding and subtract- 
ing suitable chosen pieces. The general setup for the subtraction method at NNLO has 
been discussed in [4, 5]. This setup involves three types of subtractions: 

(a) Subtraction terms between the (n + 2)-parton Born amplitudes and the (n + l)-parton 

one-loop amplitudes. These are already known from NLO calculations [6]- [10]. 

(b) Subtraction terms between the (n + l)-parton one-loop amplitudes and the n-parton 
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two-loop amplitudes. 

(c) Subtraction terms between the (n + 2)-parton Born amplitudes and the n-parton 
two-loop amplitudes. 

In this paper I discuss the subtraction terms corresponding to case (b), e.g. for one- 
loop amplitudes with one unresolved parton . I give all relevant subtraction terms for 
electron-positron annihilation into less than 4 hard partons together with their integrated 
counterparts. The restriction to e + e~ — > 3 jets and e + e~ — > 2 jets leads to simplifica- 
tions, since in these processes correlations between three hard partons are absent due to 
colour conservation [11]. Apart from that, e + e~ — > 3 jets and e + e~ — > 2 jets are from 
a phenomenological point of view the most relevant processes in electron-positron anni- 
hilation and the only ones where the necessary two- and one-loop amplitudes are known 
[12]-[19]. Valuable information for the construction of the subtraction terms for one-loop 
amplitudes with one unresolved parton is provided by the known singular behaviour of 
one- loop amplitudes in the soft and collinear limits [11], [20] -[25]. 

The remaining (and more involved) subtraction terms for double unresolved configura- 
tions will be treated elsewhere. 

This paper is organized as follows: In the next section a short summary on the used 
notation is given. Since the subtraction terms are split into several sub-pieces, this section 
is included as a "reference guide" to improve the readability of the paper. Sect. 3 reviews 
the general setup for the subtraction method at NNLO and the collinear limits of one- 
loop amplitudes. In Sect. 4 the squares and interference terms of the singular factors 
are presented. Sect. 5 reviews the soft behaviour of one-loop amplitudes and shows how 
to obtain the full subtraction terms from the known soft and collinear behaviour. The 
absence of correlations among three hard partons for processes with less than four hard 
partons is also discussed in this section. In Sect. 6 the subtraction terms are given. The 
integration over the unresolved phase space is performed in Sect. 7. Finally Sect. 8 
contains a summary and the conclusions. 

2 Notation 

In this section I give a brief summary for the labelling of the subtraction terms. The 
subtraction terms are classfied accoring to the following criteria: 

• The number of unresolved particles and the number of loops. Subtraction terms 
labelled dai'''' approximate contributions with at most n + k partons and at most I 
loops in the case where k partons become unresolved. In general they approximate 
not only matrix elements but also lower order subtraction terms. To distinguish which 
term is approximated, the notation da?^,, n is used and denotes an approximation 

of da®' k '\ The matrix element da® with I loops is identified with da®°\ This paper 
deals subtraction terms for one-loop amplitudes with one unresolved parton, e.g. it 
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provides the subtraction term 




(1) 



This subtraction term has to approximate the interference term of the one-loop am- 
plitude with the Born amplitude, as well as the integrated subtraction term da^i- 



• The particles involved in the splitting. For 1 — > 2 splittings it is sufficient to consider 
the splittings q — > qg, g — > gg and g — > qq. 

• The decomposition into partial amplitudes. The full amplitude is decomposed into 
partial amplitudes with different colour structures. The squares and interference 
terms of the full amplitudes contains therefore in the off-diagonal entries in colour 
space products of partial amplitudes with different cyclic ordering. The notation 



is used to denote such a product of two partial amplitudes with different cyclic 
ordering. 

• The decomposition into primitive parts. This follows from the decomposition of one- 
loop amplitudes into primitive amplitudes. There are three types: A leading-colour 
piece, denoted by the subscript "lc" , a piece proportional to the number of quarks Nf 
and suppressed by 1/N, denoted by "nf" and a sub-leading colour piece, suppressed 
by 1/N 2 and denoted by "sc". 

• The correlations between the hard partons. In the singular limits of one-loop am- 
plitudes with more than 4 partons there are in addition to the correlations between 
two hard partons also correlations involving three hard partons. Pieces, which only 
involve correlations between two partons (e.g. the emitter and the spectator) are 
denoted by a subscript "intr" (for "intra"). The pieces, which involve correlations 
among three hard partons are classified whether the emitter is correlated with an 
additional hard parton (denoted by "emit") or whether the spectator is correlated 
with the additional hard parton (denoted by "spec"). A similar reasoning applies 
to the singular limit of the integrated subtraction term doc^+l- Here it is useful to 
distinguish the case "intr" further, depending which two of the three partons from 
the set { emitter i, emitted particle j, spectator k } are correlated. These cases will 
be labelled "intr, i,j", "intr, j,k" and "intr,i,k". 

3 Review of basic facts and concepts 

This section first reviews the general setup of the subtraction method at NNLO and dis- 
cusses then the decomposition of one-loop amplitudes into primitive amplitudes. Through- 
out this paper I work with renormalized amplitudes. Therefore, ultraviolet renormalization 




(2) 
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is briefly discussed. Finally, the singular behaviour of one-loop primitive amplitudes in the 
collinear limit is summarized. 



3.1 General setup for the subtraction method at NNLO 

The NNLO contribution to an observable O whose LO contribution depends on n partons 
is given by 

{0) NNLO = (3) 

J O n+2 (pi,...,p n +2) da n % + J O n+1 (p 1 ,...,p n+1 ) da n 1] +1 + J O n {p u ...,p n ) da n 2) . 

The observable O has to be infrared safe. This requires that whenever a (n + r) parton 
configuration pi,...,p n+r becomes kinematically degenerate with a n parton configuration 
p'xv^p'n we m ust have 

O n+r (p 1 ,...,p n+r ) -> O n (pi,-,K)- ( 4 ) 

The subscript n for O n and do$ indicates that this contribution is integrated over a n- 
parton phase space. The various contributions da$ are given by 

j ff (o) _ (AO) * AO) \ , , 

aa n+2 — lAi+2 ^n+2J a( Pn+2, 

do? = (4» , ii? + if*il ,, a ,, '4")#, (5) 

where A$ denotes an amplitude with n external partons and I loops. d(p n is the phase 
space measure for n partons. Taken separately, each of the individual contributions da^l 2 i 
rfci+i and dcr n 2 ^ gives a divergent contribution. Only the sum of all contributions is finite. 
To render the individual contributions finite, one adds and subtracts suitable terms [4]: 

{0) N n NL0 = J \o n+2 da n %-O n+1 oda^l-O n oda^) 
+ J (0 n+l da n % + O n+1 o da { Xl - O n o da^) 
+ J (O n da® + O n o da^ + O n o da^) . (6) 

Here doi£+l is a subtraction term for single unresolved configurations of Born amplitudes. 
This term is already known from NLO calculations. The term dan' 2 "' is a subtraction 
term for double unresolved configurations. Finally, dab 1 ' 1 ^ is a subtraction term for single 
unresolved configurations involving one-loop amplitudes. The term dab 1 ' 1 ^ is the object of 
study of this article. 
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3.2 Primitive amplitudes 

The squares and the interference terms of the amplitudes in eq. (5) entangle the kine- 
matical structure with colour correlations. For a clear understanding of the factorization 
properties in soft and collinear limits it is desirable to disentangle the colour factors from 
the rest. This can be done either by introducing colour charge operators [7] or by decom- 
posing the full amplitude into partial and primitive amplitudes. In this paper I follow the 
second approach. For Born amplitudes a decomposition into partial amplitudes is suffi- 
cient. However for one-loop amplitudes it is necessary to decompose the one-loop partial 
amplitudes further into primitive amplitudes, to obtain a clear factorization in singular 
kinematical limits [24]. Primitive amplitudes can be defined as the sum of all Feynman 
diagrams with a fixed cyclic ordering of the QCD partons and a definite routing of the ex- 
ternal fermion lines through the diagram. Colour factors and coupling factors are stripped 
off from primitive amplitudes. 

Throughout this paper I use the following normalization for the colour matrices: 

Tr T a T b = -5 ab . (7) 
2 v ; 

Further, a s / (2ir) is used as the expansion parameter for the loop expansion, e.g. a n-point 
amplitude is perturbatively expanded as 

A n = c (47ra s ) n/2_1 

where c is a prefactor containig the electroweak coupling factors and C is a colour factor. 
The notation C o A n indicates that in general this contribution is a sum of terms 

Co 4° = £<f -A®. (9) 

The generic decomposition of one-loop partial amplitudes into primitive amplitudes is given 
by 

^ = ^ijc + ~Jy~ ^n,nf + JpA$ sc , (10) 

where N is the number of colours and Nf the number of quark flavours. The one-loop 
singular functions, describing the collinear or soft limit of one-loop amplitudes, follow the 
same pattern of decomposition into primitive parts: 



Sing^ 1 ' 1 ) = 




Co ^ 0) + O Co ^ 1)+ (27) 2Co ^ )+O(a 



(8) 
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3.3 Ultraviolet renormalization 

The general structure of the factorization of one-loop primitive amplitudes in singular (e.g. 
soft and collinear) limits is 

= Sing^.A^ + Sing^.Ail, (12) 

where the function Sing corresponds either to the eikonal factor (soft limit) or a splitting 
function (collinear limit). Eq. (12) holds both for bare and renormalized quantitites. 
The relation between the renormalized and the bare expression in the MS scheme for the 
singular function is 

Sing^ = SfV* Sing^ - § SingE, (13) 
where /i 2 is the renormalization scale, 

S £ = (47r) £ e- £ ^, (14) 

is the typical phase-space volume factor in D = 4 — 2e dimensions, 7^ is Euler's constant, 
and (3q is the first coeffcient of the QCD /5-function: 

A> = ^rC A - 2 -T R N f , (15) 

with the colour factors 

N 2 - 1 1 
Ca = N, Cf = -^-, T r = -. (16) 

3.4 Collinear limits 

In the collinear limit tree amplitudes factorize according to 

4°)(..., Pi ,p„...) = ^SpJit^^p,-)^-^,-)- (17) 

A 

where pi and pj are the momenta of two adjacent legs and the sum is over all polarizations. 
In the collinear limit we parametrize the momenta of the partons % and j as [7] 

, k 2 n 
Pi = zp+ kj_- 



z Ipn 

Here n is a massless four-vector and the transverse component k± satisfies 2pk± = 2nk± = 
0. The collinear limits occurs for k\ — > 0. The splitting amplitudes Split^ ' 1 -' are universal, 



7 



they depend only on the two momenta becoming collinear, and not upon the specific 
amplitude under consideration. The splitting functions Split* 0,1 ) are given by 



Splits = —uip^ip^uip), 



2 

s P nt ™ = — Hvi) ■ £ (pj) Pi ■ e (p) + £ (pj) ■ e (p) Pi ■ e (Pi) - e (Pi) ■ e (p) Pi ■ £ (pj)} » 



Sij 

Split^- = — u( Pi )f( P )u( Pj ). (19) 



One-loop primitive amplitudes factorize according to [20]- [24] 
A^(..., Pi , Pj ,...) = 

£ Split * 0,1 )( M ) ^(....p, ...) + £ Spht^fe,^) A^.^p, ...). (20) 

A A 

Here, the splitting amplitudes Split* 1 ' 1 - 1 occur as a new structure. The one-loop splitting 
amplitudes Split* 1,1 - 1 are decomposed into the primitive parts 

Split^, Splitiy>, Split! 1 /). (21) 

For the splitting q — > qg one finds [23, 24]: 

1 



splits = ^" lcr (^)"{ 

2pi -e(Pj) ; 



' 1 



-f 1 (l-z) + f : 



Splits 1 ) 



/a „a j ' «(pOiMp) » ( 22 ) 



where the coefficients /i(-z) and / 2 are 



fi(z) = ~ 2^1 ( 1, -e, 1 - e; 



and cr is given by 



e 2 \ 1 — z 

h = 2(l-e)(l-2£)' (23) 
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For the splitting g — > gg one finds 



Split 



g->gg,lc 



S-c r ( 



-? + /i(*) + /i(l-*) 



Splitf " 1 * 



2(3 - 2e) 
11 

6s 



/2 



2(Pi - Pi) • e(p) 



2e(pi) • e(pj) - 



9^99 

2pi • e(pj)2pj -e(pi) 



Split^, 



s P<iL/ 



-Si 



2(l-2e)(2-2e)(3-2e) 



2 (Pi~Pj) -e(p) 



2e(pi) • e(pj) - 



2pi -e(pj)2pj -e(pi)\ , 1 



+ — Splits 



,i) 



3e ~ r -"™9> 



Split?' 1 * 



0. 



(25) 



Finally, for the splitting g —> qq one finds 



2e(l-e)(l-pe) 



(l-2e)(3-2e) 



Splitg^gg SC 



Split^- 
2(1 -e) 
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.(o.i) 



,(l-2,)(3-2eT Split ^ + 3| Split 



(o.i) 



~ S ij 



1 e(l-pe) 



l-2e 2(1 -e) 



Split^^. 



The e-expansion of the hypergeometric function is 

2-P1 ^1, -e, 1 - e; 
The hypergeometric function can also be written as 



(26) 



(27) 



Fi 1,-£,1 -e; 



1 - z 



= r(l + e)r(l-e) 



1-2 



+ 1 - z £ 2 Fi (e, e, 1 + e; 1 - z) , 

(28) 



which is useful when studying the soft limit. 



4 Squares and interference terms 

Due to the sum over spins in eq. (17), spin correlations are retained in the collinear limit 
of squared tree amplitudes. With 

p(o,i) = J2u{p) Split^Split^ufr) for quarks, 

A, A' 
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p(o,i) = J2^(p)* Split^'Split^e^) forgluons, 



(29) 



A, A' 



and the paramterization eq. (18) one finds 



p(0,l) 



q^qg 



p(o,i) = ±_ 



g^-gg 



3 (0,1) 

g^qq 



2z 
1-z 



+ (l-p £ )(l-z) 
2z 2(1-2) 



1-2 



+ 4s(l - z) 



z 

rv A. j 



4(1 - pe)z(l - z) 



rv ft/ j 



(30) 



I introduced the parameter p, which specifies the variant of dimensional regularization: 
p = 1 for the CDR/HV schemes and p = for the FD scheme. For the interference term 
of Split (0,1) with Split (1,1) the analog of eq.(29) is given by 



pW) = J2 U (P) Split (0>1) * Split (1 ' 1} u(p) +c.c. for quarks, 

A, A' 

p(i,i) = £V(p)* SpUt (0>1) * SpUt (1>1) e"(p) +c.c. forgluons. (31) 

A, A' 

Here "c.c." denotes the complex conjugate. One finds for q —> qg 
P^lu = S,- , cr(^)"'{/ I W« + /,^[l-Ml- Z )]}-^« + 



C.C, 



p(l.l) = J_p(0,l) , 

q^qg,nf q^qg ~ u - u -> 



p(l,l) = £-l cr 

J q-^qg,sc U F 



-Si 



PK> +/2-jJ[l-Ae(l-«)]Ucc. 



7 »j 



(32) 



For the splitting g ^ gg one has 



3(1,1) 

g^gg,lc 



□ (1,1) 
g^gg,nf 



p(i,i) 
g->gg,sc 



1 



-Si 



+ / 1 (z)+/ 1 (l-z) 



3 - 2e s. 



p(o,i) 



A i 



11 

6e 

2 



/i 2 y (l-2e)(2-2e)(3-2e)sij 



_ _ p(o,i) + c c 
r g^gg ^ uu ' 



[4 - 8p^(l - z)] 



rv ft/ j 



,J_ p(0,l) , 
+ 3£ ™ ' 



0. 



(33) 
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Finally for the splitting g qq one finds 



3(1,1) 

g^qq,lc 



3(1.1) 
g-+qq,nf 



3(1,1) 
g^qq,sc 



'j 



e(l - pe) 



e 2 \l-2e 2(1 - e) 



2e(l -e)(l - pe) 
" (1 -2e)(3-2e) 



P (o,i) _ 11 p (o,i) , 
g^QQ q £ 1 q^qq 



-s-W 1 SlJ 



" 2(1 -6) (o,i) J_ ( 0) l) 

/i 2 ; 5(1 - 2e)(3 - 2e) 3^ ™ ' 



1 e(l-pe) 



l-2e 2(1 -e) 



P (0 ' 1} - +cc 



5 Soft gluons and correlations 

If a single gluon becomes soft, the partial tree amplitude factorizes according to 

A^(..., Pi , Pj ,p k , ...) = E&°>V(p i ,p j ,p k )A!® 1 (...,p i ,p k , ...), 
where the eikonal factor is given by 



2pi ■ efa) 2p k ■ e(pj) 



In the soft-gluon limit, a one-loop primitive amplitude factorizes according to 

A£\...,pi,pj,p k ,...) = Eik {0 ' 1] (pi^j^k) •A^l 1 (...,p i ,p k ,...) 

+Eik (1,1) (p i ,p J -,p fc ) • A^!_ 1 {...,p i ,p k , ...)• 

The one-loop eikonal functions are given by [24] 



EikfcV (pi,pj,p k ) 



Eikfy^ (pi,Pj,p k ) 
Eik^\p hPj ,p k ) 



(-Sij)(-S jk ) 
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Eik^fe,^-,^), 



^Eik^fe,^,^), 



0. 



(34) 



(35) 



(36) 



(37) 



(38) 



As in the tree-level case with one soft gluon, the eikonal factor is independent of the variant 
of the dimensional regularization scheme used. The square of the Born eikonal factor is 
given by 



Eik^ 1 \p i ,p j ,p k )*EikV> 1 \p i ,p j ,p k ) = 4 



(39) 
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For sub-leading colour contributions we have in addition the case of interference terms 
of two primitive amplitudes with non-identical cyclic ordering of the external partons. 
Therefore we need 

Sik , n Sn Sm 



£ Efc^W^p,) *E^\ Pl , Pj ,p k ) = 2-^ + 2-^-2 

SijSjk SijSji SkjSjl 

Eik^\ Ph , Pj , Pi yEik^\p t ,p vPk ) = 2^ + 2^-2^-2^.(40) 

SijSji ShjSjk SkjSjl ShjSji 

It should be noted that the product of the two eikonal factors, which share only one or no 
hard parton decomposes into squares of eikonal factors. From eq. (39) and eq. (40) one 
obtains at NLO all subtraction terms by partial fractioning 

S S ' + (41) 



SijSjk Sij^Sij + Sjk) Sjk(.Sij + Sjk) 

and by extending then each term on the r.h.s of eq. (41) to the full dipole splitting function 
[26]. For one-loop amplitudes with one unresolved parton, additional e-powers of eikonal 
factors occur. For example 

/i 2 (-Sifc) 



E V^\ Pl , PvPl )*E^\ Pi , P j, Pk ) 



2 



/j 2 (-gjfc) 

~ s ij)(~ s jk) 



$ik $il Sik Ski 



_^ Sil Ski 



SijySij + Sjk) Siji^Sij + Sji) SjkyS^j + Sjk) Sjk\Sjk Sji) 

(42) 



Sj^Sij + Sji) SjiiSjk + Sji) 

Here one can distinguish three types, depending on the hard partons which occur in the 
£-power of the eikonal factor. Contributions of the form 

V 2 (-Sik) Y s ik ^ 



s ij)( Sjk) J Sij(Sij + Sjk) 

involve in addition to the emitted soft parton j only the emitter i and the spectator k. 
These contributions are labelled "intr" . Terms of the form 



( s ij){ Sji) J Sij(Sij + Sjk) 



(44) 



are labelled "emit", since the eikonal factor which occurs to the e-power corresponds to 
the emission of a soft parton j from an antenna formed by the emitter % and an additional 
hard parton I. Finally, terms of the form 

Z 1 {—Ski, i - ih- 



( s kj)( Sji) J Sij(Sij + Sjk) 
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are labelled "spec" , since the eikonal factor which occurs to the e-power corresponds to the 
emission of a soft parton j from an antenna formed by the spectator k and an additional 
hard parton I. 

In general, these correlations with an additional hard parton / have to be taken into 
account. However, for the processes e + e" — > 2 jets and e + e~ — > 3 jets these correlations 
are absent. This is obvious for the NNLO calculation e + e~ — > 2 jets, since there are only 
two hard partons. It also holds for the NNLO calculation e + e~ — > 3 jets due to colour 
conservation. The vanishing of the accompanying colour factor is most easily seen when 
working within the formalism of colour-charge operators. In the formalism of colour charge 
operators it is convenient to introduce for an amplitude with n + 1 partons an orthogonal 
basis of unit vectors \c±, c 2 , c n+ i) in the colour space. An abstract vector \M. n+1 ) in this 
colour space is then defined as the projection of the amplitude A n+ \ onto this basis: 



At+i = (ci,c 2 , ...,c n+1 \M n+1 ) 
The vector \Ai n +i) is expanded in the coupling a s as 



(46) 



\M 



n+ll 



= c (to,) 1 -""^-!- 1 '" 



*&,) + (" 



l n+l 



V27T 



M%) + 0(ai) . (47) 



In the soft limit where the momentum of parton j becomes soft, the n+l-parton amplitudes 
behave as 



M 



(o) 

n+l 



(or 



A<Si) = V^S S~^( Pj ) [j(°>(p,) |M«> + (g) jW(p,) |M°)>] , (48) 
where the soft currents are given by [11] 



j Pj 

J^ (1) fe) = 



r(i + £ )r(i-5)-^r bc ^T^ 



H y—Sik) 

~ s ij) (~ s jk) 



Pi 



Pk 



Pi ■ Pj Pk ■ Pj 



_^J-(0)( p ). 

2e vw 



(49) 



The colour charge operators Tj for a quark, gluon and antiquark in the final state are 

quark: (...g i ...|7?-|...q J -...), 



gluon: (...g c ...\if 



cab I 



■5 



antiquark : (-^) l-ty--)- 



Colour conservation reads 



E T *K 0) > = °- 



(50) 
(51) 
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From eq. (49) it follows that terms labelled "spec" or "emit" are accompagnied by a colour 
factor 

if abc T^T b k Tl (52) 

This colour factor equals zero for processes with less than four hard partons. Therefore for 
the major applications e + e~ — > 3 jets and e + e~ — > 2 jets these correlations are absent. 
Similar considerations apply to the approximation dajj'jj n of the integrated subtraction 

term dot^l- The integrated subtraction term da^^l can be written as 



aa n+l — g 7r 2 \ JVl n+l 



M%). (53) 

with the insertion operator 



j_ v t?t- sjH^f ( Slm y 

I--^^p^ m , Z lm _ 1 — 1 V,. (54) 

The singular factors V; are given for a quark by 

V q = CpVjM,, (55) 

and for a gluon by 

= C A V^ 9 + N f V^% (56) 

The functions V*- ' 1 -* will be given in eq. (95). In the soft limit the amplitude Ai^iJ 

factorizes as in eq. (48). If we denote the emitter by i, the emitted soft parton by j and 
the spectator by k, we can distinguish the following cases: 

• l,m e {i,j, k}. In this case we obtain terms of the form 

\ — £ / \ — £ / 



These cases are labelled "intr, i,f\ "intr, j, h" and "intr, i, fc", respectively. 



(57) 



Either 1 or m in A;}, and the other one distinct. In this case one obtains terms 
of the form 

Sii \ ~ £ s ik ( s kl \~ £ s ik ( Sji\~ £ s ik ^ 

Z^ 2 / Sij(,Sij 4" Sjk) \ Z^ 2 / ^ij{sij 4~ Sjfc) y/i 2 y Sy(sjj + Sj&) 

With an argument similar to the one used for the soft behaviour of one-loop ampli- 
tudes, one can show that the colour factor for the case sjf vanishes for amplitudes 
with less than four hard partons. Therefore only the first two cases have to be 
considered. They will be labelled "emit" and "spec". 

l,m $l k}. This case requires at least four hard partons and is therefore not 
relevant to e + e" — > 3 jets and e + e~ — > 2 jets. 
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6 Subtraction terms 

The approximation terms dctn' 1 ^ and doi^+l are written as a sum over dipoles [4, 7]: 



(59) 



topologies T 



topologies T 



The splitting topologies are here just 1 — > 2 splittings. A dipol is constructed from ampu- 
tated amplitudes, evaluated with mapped momenta and dipole splitting functions. By an 
amputated amplitude I mean an amplitude where a polarization vector (or external spinor 
for fermions) has been removed. I denote amputated amplitudes by \A). If the amputated 
parton is a gluon, one has 



\A n (...,p,...)) = 



d 



A n (...,p, ...). 



If the amputated parton is a quark, one has 

\A n (...,p,...)) = 



d 

du{p) 



A n (...,p, ...). 



(60) 



(61) 



The dipole subtraction terms are obtained by sandwiching a dipole splitting function in 
between these amputated amplitudes: 



VS$F) = |c| 5, (47ra.) B C y (^ 1;l (...,p e , ...) V^(T) ^ 1;j (...,p e , ...) 



|(0) 



(o,i), 



i(0) 



.(62) 



Here, c denotes an overall prefactor and C is a colour factor. The subtraction term dan 
approximates 



(o,i) 



E 



|q (47rct s ) 
8tt 2 



4(0) ' 4W , a\l) ■ a\u) , 4W 7-. . 4 



l(l) 



|(0) 



>(0) 



,(0) 



# n+ i.(63) 



Xij is obtained from integrating da^+l over the unresolved phase space. The subtraction 
terms for this contribution are written as 



P (i,i) (T) 



I I 2 I A V 

\c\ (47ra s ) 



8tt 2 



p(o>i)(7 1 ) 



1(1) 



p(o,i)(y) 



4 (0) 



(0) 



1(0) 



4 (o) 

Alii 



4 (0) 



(64) 



Here, in the first and second term the NLO splitting function T^ 0,1 ) is sandwiched between 
A ( n ] and A { n\ If we recall the factorization formula eq. (12) for one-loop amplitudes 



4(1) 



SingW.^ + SingW-Af, 



(65) 
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then the first two terms of eq. (64) approximate exactly the singular part proportional to 
Sing^ - 1 . The third and fourth term in eq. (64), which involve the splitting function Vj{'y, 

approximate the singular part proportional to Sing*- 1 -*. Finally, V^'^ in the last term of 

eq. (64) is an approximation to da^+i- 

The amputated amplitudes are evaluated with mapped momenta. The mapping of mo- 
menta relates a n + 1 parton configuration to a n parton configuration (or a n + 2 parton 
configuration to a n + 1 parton configuration). Such a mapping has to satisfy momen- 
tum conservation and the on-mass-shell conditions. Furthermore it must have the right 
behaviour in the singular limits. Several choices for such a mapping exist [7, 27]. One 
possible choice relates three parton momenta of the n + 1 parton configuration to two 
parton momenta of the n parton configuration, while leaving the remaining n — 2 parton 
momenta unchanged. This mapping is given by 



Pe = Pi+Pj 



Pk, 



Ps = 



i-y 



Pk- 



(66) 



Pi, pj and pk are the momenta of the n + 1 parton configuration, and p e and p s are the 
resulting momenta of the of the n parton configuration. The variables y and z are given 
by 



$ik 



y = , z = 

Sijk s ik ~r Sjfe 

It is also useful to introduce the vector k±, defined by 



(67) 



k± = (1 - z)pi - zpj - (1 - 2z) 



Pk- 



The dipole splitting functions V^ 0,1 ^ read 



(68) 



T>(0,1) 



-p(o,i) 

9^99 



? (0,1) 

9^m 



2 ^ 
P 

$ijk y 

2 1 

&ijk y 

2 1 

Stjk y 



Zl 



-9^ 



y) 

2 



:9 



1 - z{\ 
4r(l - 



2 + (l-pe)(l-y)(l-z) 



-2)+(i-P£) 4r(1 y)2 y\k\ 



-y) 
y? 



ySijk 



yS-ijk 



(69) 



For the q — » qg splitting, the quark polarization matrix j> e is included in the definition 
of the dipole splitting function. For the g — » gg and g — > qq splittings, the full splitting 
function is shared between the dipoles with spectators k and k', where k and k' are the 
partons adjacent to the pair in the colour ordered amplitude A^°\...,k',i,j,k,...). 

For the g — > gg splitting the dipole splitting function is chosen in such a way that the 
singularity when parton j becomes soft resides in V g ^ gg (pi,Pj;pk), where as the singularity 
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when parton i becomes soft resides in V g ^ gg (pj,Pi;pk>)- The parameter r is either chosen 

as 

1 

r = 2' 

or as 

2pk>Pe 



(70) 



r = 



2p k ,p e + 2p e p s ' 



(71) 



where k' is the other parton adjacent to the emitter % in one of the amplitudes. The 
g — > gg and g — > gg splittings also involve spin correlation through the spin correlation 
tensor k^k 1 ^. Note that the spin correlation tensor /c^/c^ is orthogonal to p e and p s : 



2p e /c_L = 2p s k_L = 0. 



(72) 



Further note that the contraction of p e into an amputated amplitude vanishes due to gauge 
invariance. 

The parameter p specifies the variant of dimensional regularization: p — 1 corresponds to 
the CDR/HV schemes and p = to a four-dimensional scheme. 
The splitting functions f° r the splitting q qg are given by 



"T> (1,1) = 

' (1,0) q^qg,lc,corr °e c r 
X <gi,corr(y,z)V^' 1) 



i 



p. 



(0,1) 



(1,0) q->qg,nf g £ g-+g<?' 

-p(l,l) ' _S «^ 

7 (1,0) g^?g,sc °£ c r 



X 



- 9l,intr(V, 1 - Z) 



2 1 



^S + /2 i> e [l-pe{l-y){l-z) 



Sijk y 



The label u corr" takes in general the values {intr, emit, spec}. However for e + e 
or e + e~ — > 2 jets only pi^,. is relevant. The function gii ntr is given by 



(73) 
3 jets 



gi,intr(y,z) 
1 

~72 



r(i + e )r(i - £ ) (— + 1 - (i - y y z z 2 f x ( e , e , 1 + e; (1 - y )(i - *)) 



(74) 



For the numerical integration over the (n+ l)-parton phase space we need the ^-expansion 
of the hypergeometric function (e, e,l + e;x): 



2 F 1 (e,e,l + e;x) = 1 + e 2 Li 2 (x) + 0(e 3 ). 



(75) 
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For the splitting g — > gg one has 



T>(1,1) 

(1,0) g-^gg,lc,corr 



1 



~$ijk 



+ 9i ,corr 

(y, z) + gi yin tr{y, 1 - -2) 



T>(0,1) 



9~*99 



3-2 £ /2 4 fe y2 



T 3 , 



(1,0) g^gg,nf 



X 



2 r 

s lk y 2 



2 r 



— e 

~°ijk \ — F 



2(1 " ^) 



8p£(l-y) 5 



11 



2(1-2) 



(l-2e)(2-2e)(3-2e) 
8pe(l-y) 2 



i 

, ±_ <p(0,l) 



(1,1) 

(1,0) g^gg,sc 



0. 



(76) 



Again, only the label corr = intr is relevant for the applications e + e — > 3 jets and 
e + e~ — > 2 jets. Finally for the splitting g — > gg one finds 



-p(i.i) _ c-i„ 

' (1,0) fl -M,ic - D e C r 



II 2 ' " 



+ 9i,intr(y, z) + gi,i„tr(y, i - 2) 



1 f 1 e{l-pe) 2e(l-e)(l-pe) 



e 2 Vl-2e 2(1 - e) (l-2e)(3-2e) 



-p(o,i)_ _ n -pi 



11 



>(0,1) 



T>( 1 ' 1 ) -_Q-1„ 
' (1,0) g^qq,nf ~ D e C r 



2(1 -e) 



') y e(l-2e)(3-2e) 



-p(o,i)_ 1 _ -p 



(0,1) 



V, 



(1,1) 

(1,0) g^qq,sc 



~$ijk \ 



2 - 1 y ~ 



1 e{l-pe) 



l-2e 2(1 -e) 



P. 



3e 

(0,1) 
g^qq- 



q^qqi 



(77) 



Let us now turn to the subtraction term da^'^ n and the splitting functions P^'j'. The 
splitting functions P^'ij are gi ven by 



(1,1) 



5 , £ 1 (4tt) £ /ggA £ ^l- e V (0,1) P (0,1) 



(0,l)a->6c ~~ r(l — £) \ fl 2 ' L"'-<"-'A.7- -)\ y u'—b'c'>u^bc- (^ 

where a' — > fe'c' is the splitting, which has been integrated out in the subtraction term 
da^+l- The label "corr" takes the values u intr i,j", u intr j,k" , il intr i,k", "emit" or 
" spec" . If r = 1/2 is chosen for the subtraction terms da^l, then V$^ b , c , is independent 
of any kinematical variable. The function h is given by 

hintr, 1 Z) 2/ 5 

h in tr,j,k(y^ Z ) = ( 1 -y)( 1 -2), 

h in tr,i,k(y^ z ) = (1-^)2, 

2p e p« 



h e mit(y,z) = (l-y)z- 



$ijk 



h special Z) 



Sijk 



(79) 
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Here, p\ is the additional hard momentum appearing in eq. (58). 



7 Integration of the subtraction terms 

The phase space measure for n + 1 particles in D dimensions factorizes according to 

J d(fi n +i(P > Pi, pj, pk, ...) = J d<f> n (P ► ...,p e ,p s , ■■■)d(f) unreso i ve( i{sij, Sjk, Sjfc, Q^_ 2 ). 

(80) 

where f^_ 2 denotes the dependency of the unresolved phase space on the D — 3 remaining 
angles of particle % in D dimensions. In the absence of correlations between three hard 
partons this dependence is trivial and the unresolved phase space is given by 



dcf). 



'unresolved 



(4tt) £ - 2 

The variables y and z are given by 



s ijk f- £ J dyy-%1 - y) 1 ^ j dzz'%1 - z)' £ . (81) 



V=^, z=—^. (82) 

The collinear limit is obtained by y — > and in the soft limit pj — > one has in addition 
z — > 1. The splitting functions for g ^ gg and g ^ qq involve the spin correlation tensor 
k^k 1 ^. The integral over the spin correlation tensor can be written as 

c*M + c 22P K + C23 tor + rtr® + <W". (83) 

Using 2k±p e = 0, 2fcjj> s = 0, pi = and p 2 s = this reduces to 

d<f> unr esolve d f{y,z)k»kl = ~C 24 (-g^ + 2 ^±^- ). (84) 

Due to gauge invariance only the term — C*24(— g^ v ) will give a non-vanishing contribution. 
C24 is obtained by contracting with g uu , 

If n v 

C24 — 7^- -j^pj J d(p unreso i vec [ f{y, z) gavk^k^-, (85) 

and therefore 

d4>unresolved f(Vi z ) k^_kj_ = 

' S " A ' / g ^ + 2 P " P l + P ° K ) [d<j> unresolved f(y,z)yz(l-z). (86) 



2(1 -pe)\ ■ 2 PeP , 
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The integration over the unresolved phase space can be reduced to three basic integrals. 
The first one, 

1 1 



J dy y a (l - y) l+c+d J dz z c (l - z) d [1 - z(l - y)f 



r(l + a)r(l + c)r(l + d)T(2 + a + d + e) 

r(2 + a + d)r(3 + a + c + d + e) ^ ' 

is sufficient for the integration of the dipole factors V^ 0,1 \ Note that the integral in eq. 
(87) depends only on four free parameters a, c, d and e. The exponent of the factor 
(1 — y) is fixed by these parameters. At the one- loop level there are in addition integrals 
over the functions g±(y,z) and gi(y, 1 — z), which in turn involve an integration over a 
hypergeometric function. These are given by 
l l 

J dy y a (l - y) 1+c+d J dz z c (l - z) d [1 - z(l - y)] e 2 F, (e, e; 1 + e; (1 - y)(l - z)) = 
o o 

r(i + a)r(i + c )r(i + e) 



T(e)T(e) 



oo 



, T(j + ejTjj + ejTjj + 1 + d)T(j + 2 + a + d + e) 

F(j + l)r(j + 1 + e)T(j + 2 + a + d)T(j + Z + a + c + d + eY 



and 

i i 

J dy y a (l - y) 1+c+d J dz z\\ - z) d [1 - z(l - y)f 2 F X (e, e; 1 + e; (1 - y)z) = 

o o 

T(l + a)r(l + d)T{2 + a + d + e)T{l + e) 

T(2 + a + d)T(e)T(e) 

x v r(j + e )r(j + e )r(j + i + c) 
^r0' + i)r0' + i + e)r0' + 3 + a + c + d + e)' 1 ' 

The results are proportional to hypergeometric functions 4F3 and 3F2, respectively, with 
unit argument. These sums are then expanded into a Laurent series in e with the techniques 
described in [28, 29]. For the symbolic calculations the computer algebra system GiNaC is 
used [30]. To simplify the resulting expressions, the following formulae are useful: 

2 4 4 4 4 

_ 7T Z _ 7T 7T 7T _ 7T 

G = ¥' G= 90' Cl > 2 = G ' Cl ' 3= 3^0' C2 ' 2= 120' Cm ' 2 = 90- (90) 
For the multiple zeta-values the notation 



Cm k ,...,m 2 ,mi / . -mi -mi -rrik 

(91) 



il>i2>--->ifc>0 



is used. 
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7.1 Integration of the subtraction terms for single unresolved 
contributions 

We write the integrated dipole factors as follows: 



8n 2 S. 



v/. 



dffrunresolved T^ a ^l) C {Pii Pj> Pk) 



S £ (4tt) £ ( s ijk \ ( o,i) 

2 T KUte + S au S e terms, 



(92) 



r(l-e) \fi 



where T is the appropriate polarization sum tensor, e.g. for quarks one has 



(93) 



whereas for gluons one has 



T = 'E(eZ(k,n)yeKk,n) = - ! r + 



k ■ n 



n 



(k • n)' 



(94) 



In the case for an emitting gluon, there are additional terms indicated by "gauge terms" , 
which vanish when contracted into an amplitude. Integration of the subtraction terms for 
single unresolved contributions yields: 



y(o,i) 
q^qg 



y(o,i) 

9^99 



V, 



(0,1) 



9^11 



r(-e)r(-2e)r(i- e ) 

r(l - 2e)r(l - 3e) 
17 1 



Ae 



l-3e 



2(1 -A*) 



l-3e)(2-3e) 



+ 



585 45 17 2 1 2 



99 7 3 2 



11 



+ ^rP- ^ - -*P^ - 12 Cs - ^ttvt 4 e 2 



16 8 



8 



r(-e)r(-2e)r(i- e ) 

r(l - 2e)r(l - 3e) 



2 + 



4e 



120 

4 e(l-e) 



+ 0(e 3 ), 



-r 



l-3e 3 (l-3e)(2-3e) 



7 1 



15 



/ 93 9 7 2 8 > 11 .. 

+ (54 - -r - 3n 2 + -r^ - 16 C s + -r( 3 - ) e 2 

+0(e 3 ), 
T(-eni-e) 2 



T(2 - 3e) 



1 2 



2 3(l-/oe)(2-3e). 



1 3 7 



9 15 



1 2 1 2 



o + q r J o + fi r + 77 + [ -77 + — r + -pr + -n - -m ) e 



27 93 21 3 2 
+0(e 3 ). 



2 ' 4 '2' 



(95) 
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Since p takes only the values and 1, p 2 = p was used to simplify the ^-expansion for 
v (o,D_ 

y g^>qq- 

Note that these results are sandwiched in eq. (64) between a one-loop amplitude and 
a Born amplitude and are therefore needed to 0(e 3 ). 



7.2 Integration of the subtraction terms for contributions with 
one virtual and one real unresolved parton 

We first consider V^q). We write the integrated dipole factors as follows: 



8-tt 2 q- 1 !, 26 / rlrh -pCM) 
o/l <J £ /J, I Uj^unresolved ' (i q) 

For the splitting q — > qg one finds 

,(i,D _^J_ 31 

4e 4 



^ ) T V (io) + S au S e terms - ( 96 ) 



V, 



(1,0) q—>qg,lc,intr 
\ 



12e 



3 



r(i- £ ) 2 V a* 



51 1 5 2 11 \ 1 



145 9 15, 11 r 



72* + T C3_ 4 



12 e 



4' 12 
1663 233 



6 



151 55 
-6--24 P 



16 



107 2 5 2 356 , 
24 P+ l6- 7r + 12 P7r + ^ C3 



1 4 187 r 11 r 55 2r ll r2 11 r3 . 



1 3/29 



1 vr 2 \ 1 



139 11 o 15, 

p + 7T H C3 

8 8 F 2 s 



(1,1) 



+ 0(e), 
1 5 



1 I 1 \ 1 /17 1 5 2 1 _ 1 _ 2 \ 1 
(i,o)™n/ = ^+U+3 L )^+(7i + 6P-^ +2 L+ 6 L 



33 7 5 2 23, 17 T 1 5 2r 1 .2 1 r 3 ^, x 

+ — + — p 7T 2 Cs + — ^ + ~pL n L + -L + — L + O(e), 

8 12 H 24 9 s 12 6 H 36 4 18 V ; ' 

11) k 2 \ 1 /35 3 7T 2 „ , \ 1 or 15 35 2 1 9 



24 



12' 



21 



17 



5 7T 2 \ 1 35 3 7T 2 



8 6 



I + T + 8^-T- 7C3 



(97) 



Here, L = ln(sjj fc /// 2 ). For the splitting g ^ gg one finds 



(1,1) 

(1,0) g^gg,lc,intr 



1 / 17 1 \ 1 

+ 77" +7f U 



45 4 



6 6 



29 16 7 2 11 r \ 1 
-77- + ^ + t^tt - -^L — 



3 



9 



12 



6 



949 199 2 1 2 29, 11 11 11 r2 \ 1 26875 

+ -43 H r H 7T 2 vr 2 r H (3 L H rL L 193 H r 

108 72 3 2 s3 3 18 12 J e 648 



7T 2 - 



I 415 
H or H 

18 P 36 

II r2 11,3 • 
+ 7777?"£ - 7777^ + W 



499 
216 



36 



36 



HI 7 , 7 4 n - r 77 r 55 2 r 11 r2 

*'r + ^3 - S( 3 r + jg** - 11L + -rL + -^L - — L 2 

1 / 1\1 / 7 tt 2 \ 1 359 
— + -1 + -r — + -6 + -r + — - - 32 + — r 



Ae 3 



6 



6 



54 



22 



5 2 5 2 29/ 



(1,1) 



1 



(1,0) g->gg,nf 3^3 



+ 



+ 0(e), 



9 



+ 2 - -r 

9 



36 



-7T 



-L 



9 



-rL 



Mi 

6 / £ 



1 1 5 2 5 2 

+6 H r H — pr 7r H 7r r 

36 9^ 18 108 

+ ^L 3 + mir + 0(e), 



(1,1) 



0. 



(1,0) g^gg,sc 
For the splitting g — > gg one finds 



^ Cs + 2L - -rL 
9 s 18 



36 3 18 



(98) 



V, 



(1,1) 

(1,0) g^qq,lc 

rL 
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1 


/247 


1 121 


(I 


r 

36 J 


- + 

e 2 


(n 4 


" 24 P ~ ~54 


: + 


3581 


23 


7339 


329 




72 P " 


r 


pr — 

108 F 


216 


648 



7T 



r — 
119 



—or 1 — 7r r 

9^ 6 9 



11 

12' 



L 



18 
14 

+ ^7(3^ + ~r L ~ T^rL - —prL + — L z - —rL z + in 



2 139 2 1 2 > 

7r H 7r r H — 7r or — 7C 3 

144 216 9 1 S 



11 



77 



11 



11 



11 



4 
11 



36 
1 



18 



24 



36 



1 Ml 49 

I r - H 

24 12 e 72 



7T 



H p r H or 1 — 7r r 

24^ 108 18^ 6 9 



+ 0(e), 



V, 



(1,1) 

(1,0) g^qq,nf 



4 7 11 

r L H — rL 

9 18 6 9 



151 895 7 11 2 11 2 

H r pr 7r H 7r r 

54 324 18^ 72 108 



1 r 7 r 1 r 1 r2 1 r2 • 

— L H rL + -prL L H rL 2 + ?vr 

2 18 9 r 12 18 



1 1 

6 ~ 9 r 



1 17 7 1 

— I r pr 

e 18 9 9^ 



+ 0(e) 



(i,i) 

(1,0) g^qq,sc 
25 



1 1 

1 — r 

4 6 



+ 



11 41 1 

1 r H — pr 

8 36 Q H 



51 



71 



H pr H 7T r 

18^ 3 9 1 



1 219 



5 



41 



641 
108' 



4441 212 11 2 41 2 2 2 

r H pr H 7r 7r r 7r pr 

162 27 1 6 27 9 r 



+5C3 - y C3^ + «tt 

25 7T 2 7T 2 ' 

H pr H r 

18^ 4 6 



4 ' 6' 
i 0(e). 



11 41 



— - + -r 1 — + — — + — r + -pr -p + — — r 



8 36 



1 51 



641 



108 



(99) 



Again, since p only takes the values and 1, p 2 = p was used to simplify the results. 
Let us now consider the quantities V(oi), obtained from the integration of the sub- 



traction term da^'^. The corresponding splitting functions V^'^ were given in eq. (78) 
by 



(i,i) 

(0,1) 



,(i,i) 

(0,1) a^bc 



r(l-e) \p 



(100) 
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If r = 1/2 is chosen for the subtraction terms da^+i, then V^^ b , c , is independent of any 
kinematical variable. In this case we define the kinematical functions V/q'^ by 

8n 2 S; V 2£ / d<t> unresolved = "* T vjg VS C , + gauge terms, 

(101) 

e.g. V^'^ b , c , is factored out. For two-parton correlations we have to distinguish the cases 
u intr,i,j", u intr,j,k" and u intr,i,k". One finds 

,,(ii) 1 3 11 1 5 2 /21 5 5 2 „\ 

VfcVj = ^ + + T + 4 P " 18" + U + 4 P " 8" " 6C V ' 

/ 11 55 „ 5 9 27 A 113 A 9 ,^ „. 
+ 41 + — p 7r 2 pvr 2 Cs 7T e + (e 3 , 

,,(H) 1 fi 1 \ 1 a 8 5 2 fir 217 5 2 5 2 

V W' - i = ^ + { 1 -r)-e +A -9 r --8 n+ { 16 -lA r - r + 36™ 
-6Ca) e + (6 4 - - yvr 2 + grvr 2 - 18C 3 + 3r( 3 - e 2 + 0(e% 

,,(H) f 1 M 1 n 8 1 f „ 217 19 5 2 

V (0 ;i) = (^-4 + ^ ) - - 1 + + ^ + l^ 4 + - 54 r + I8^ r + 24" 

5 2 \ f ^ 1379 137 5 2 20 2 5 2 9 ^ 0^2 
-36 rn ) " + I' 16 + ^T r + Yl pr + 6 n ~ 27™ " 36^" + 2 Cs " 3<3 ) £ 

+0{e 3 ). (102) 

Vg^™^ = i~ 2 + I + 6 + \ P ~ \? + ( 24 + 2 P - \n 2 - 12(a) e 

+ ( 9 6 + 8p - 5rr 2 - l p7 r 2 - 27C 3 - J^vr 4 ) e 2 + 0{e\ 

,,di) 2 / rt 1 \ 1 n 13 5 , / 160 5 2 5 2 
VL'V . , . , = — + 2 - -r - + 8 r - -7T 2 + 32 r - -tt 2 + — m 2 

(0,1) g^gg,%nbrj,k 3e 2 \ 3 J E 9 9 \ 27 3 18 

/ 1936 20 2 65 2 , , H3 A 2 <u 

-12C 3 ) e + f 128 - —r - -n 2 + -rvr 2 - 36C 3 + 6rC 3 - — ttM e 2 + 0( £ 3 ), 

,,di) / 1 1 A 1 13 1 / 160 16 5 2 
VL'V - • , • . = h -r 2 H r + -pr + -8 H r H or H 7T 

(0,1) g^qq,intr,3,k \ 2 3 J £ 9 3^ \ 27 9 H 12 

5 A ^ / Q9 ^ 1 93 6 2 08 5 2 65 2 5 2 , q a ,,) 2 
--rTr j e+ (-32 + — r+ — pr + -n --rvr - -prvr +^-^6 

+C(e 3 ). (103) 

,,di) 1 3 11 1 5 2 / 5 5 2 > \ 

TO™,**** = ~ 2 + ^ + T + 2 P " 6* + ( 21 + 2 P ~ r 18(3 ' 
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^S^,, 4 + - r) \ + 8 - t r - + ( 32 - ™ r - + s- 2 
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1936 20 2 65 2 „ , 113 A 2 ^, o x 

-18C 3 ) e + f 128 - — r - — tt 2 + -ryr 2 - 36C 3 + 6rC 3 - — ttM e 2 + 0(e 3 ), 
/ 1 1 \ 1 13 1 / 160 16 5 2 

V (0,l) ™intr,i,* = [~2 + 3 r J i - 2 + Y r + 3^ + ( v " 8 + ^ + + 12* 

5 2 \ ( nn 1936 208 5 2 65 2 5 2 rt> „ A \ 2 
-^- 2 J e + (-32 + —r + — pr + -vr 2 - -rvr 2 - -prvr 2 + 9 C 3 - 6r C s ) 

+(D(e 3 ). (104) 

For the NNLO calculation for e + e _ — > 3 jets it remains to treat the cases V(o'i) emii and 
"^(o ij spec These are relatively simple and given by 

V (l,l) = ( 2 PePl\ 

(0,1) a^bc,emit \ j (0,1) a^bc,intr,i,k'> 

v (i,i) _ f^PsPiY 5 v (o,d fl05 x 

If the variable r in V^?!^, is chosen according to eq. (71), then in general V„?^ 6 , c , will 
depend on kinematical invariants. This implies that V^^ b , c , cannot be factored out from 
the integration over the unresolved phase space, as it was done in eq. (101). The results 
for this case are not listed here, but can be obtained easily. 



8 Summary and conclusions 

In this paper I considered the subtraction terms for one-loop amplitudes in the case where 
two partons become collinear or where one parton becomes soft. These subtraction terms 
are relevant for precision calculations at next-to-next-to-leading order. For two- and three- 
jet production in electron-positron annihilation all necessary subtraction terms have been 
given and the subtraction terms have been integrated analytically over the unresolved phase 
space. Once the corresponding subtraction terms for double unresolved configurations have 
been worked out, the method for the cancellation of infrared singularities is complete and 
it should be possible to extend the exisiting numerical programs for NLO predictions on 
e + e _ — > 4 jets [31] -[34] towards NNLO predictions for e + e~ — > 3 jets. 
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